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Abstract 

We study the quantization of three-dimensional many-body systems in rotating coordinate 
frames defined implicitly by frame conditions. We carry out the elimination of orientational 
degrees of freedom in general, giving the Hamiltonian for the ./V-particle system in a broad class 
of body frames in terms of frame conditions and internal coordinates. We obtain several forms 
for the kinetic energy operator and compare them to related expressions in the literature. 

1 Introduction 

The problem of separating the dynamics of quantum many-body systems into collective rotations 
and internal motions leads to their quantization in rotating frames. We refer to a rotating frame as 
a body frame when the components of the total angular momentum operator in a space-fixed frame 
and in the rotating frame satisfy the same commutator algebra as in the rigid body problem. In the 
latter case the body frame is essentially unique, up to time-independent rotations and symmetry 
transformations of the rigid body. For a general iV-particle system there is a large freedom to 
choose a body frame. It is thus of interest to study the quantization of many body systems in a 
class of body frames as wide as possible. 

In this paper we study the quantization of three-dimensional many-body systems in rotating 
coordinate frames defined in implicit form by frame, or gauge, conditions. We carry out the 
elimination of orientational degrees of freedom in general, giving the Hamiltonian for the A^-particle 
system in a broad class of body frames in terms of frame conditions and internal coordinates. We 
obtain several forms for the kinetic energy operator and compare them to related expressions in the 
literature, showing how the coefficient functions are fixed by the frame convention, through frame 
conditions, and internal coordinates. In the case of linear frames and body-frame coordinates our 
results reduce to those previously obtained in 

The generic Hamiltonians discussed here can be applied to specific physical systems by choosing 
internal coordinates appropriate to the system under consideration. Those physical problems 
include, for instance, the determination of molecular rotation-vibration energy levels and their 
wave functions and scattering problems in molecular, atomic and nuclear physics. There is a 
vast body of literature on the quantization of many-body systems which we do not try to summarize 
here. Broad reviews relevant to the point of view adopted in this paper are given in pQQ]. 

* abouzas@mda.cinvestav.mx 
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The outline of the paper is the following. In the next section we discuss several technical 
issues related to body frames and the frame conditions defining them that are needed in order 
to obtain the Hamiltonian and quantum inner product for a many body system. Those include 
the form of admissible frame conditions and their reparameterizations, the body-frame angular 
momentum, and internal coordinates. In section|3]we derive the body-frame kinetic energy operator 
in standard order in terms of internal coordinates and frame conditions. The form of wave functions 
referred to the body frame and their inner product is discussed in section^ where we also give the 
representations of the total angular momentum and kinetic energy operators both as irreducible 
matrices and in rigid rotator form. In section we give two alternate forms for the Hamiltonian 
and discuss their equivalence with the standard ordered form given in section |3J Also, we make 
contact with the gauge-field formalism of [Hj by locally expressing gauge fields in terms of frame 
conditions. Examples with N = 3 and 4 are briefly examined in sections El and [7| as verifications 
of the formalism of the previous sections, and the results compared with those from the molecular 
literature. In section El we give some final remarks. 



2 Preliminaries 

We consider rotating frames whose definition depends only on the coordinates of the particles and 
not on their velocities, nor on the angular velocity of the frame itself. Those frames can always be 
defined implicitly by imposing conditions on the position vectors of the particles. Most of the body 
frames commonly used in the literature belong to this class, as illustrated in sections El and below 
with two familiar examples. Another well-known example is the Eckart frame We do not 

impose any restriction on the form of frame conditions, provided they fix the frame uniquely. The 
particular case of frame conditions depending linearly on the particles coordinates was considered 
in from the point of view of gauge invariance. 



2.1 Frame conditions and their reparameterizations 

The frame conditions defining the body frame must fix its six degrees of freedom. The translational 
degrees of freedom are fixed by choosing the center of mass frame. Thus, the frame conditions take 
the form, 1 

N 

C{{r}) = J2m a r a = , G a {{r}) = , a = 1,2, 3, (1) 

a=l 

where Q a are three conditions fixing the orientational degrees of freedom. We denote by {r} = 
{r%, . .. ,rjv} a generic configuration 2 of the iV-particle system. Since lab- frame configurations are 
not restricted, any {r} can be a lab-frame configuration. Those configurations satisfying the frame 
conditions are denoted by {R} . Thus, C({R}) = = Q a ({R}) and such a set {R} of N position 
vectors can be a body-frame configuration. We introduce also the following notations, 

|£ = *,«*» = Ej>iU, *1 (WH = ,2) 

a=l J a=l 
1 The letters a,b,c,d are used for non-tensorial indices, as in Q a . Summation over those indices and their ranges 

of variation are always explicitly indicated. We only use the summation convention for tensor indices, denoted by 

latin letters i, j, k,l, . . ., which always run from 1 to 3. Greek indices number particles. 

2 A configuration of the system is actually {r,r}. We refer to {r} as a configuration here for convenience. 
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which will be used throughout the paper. 

In order for the conditions Q a = to fix the orientational freedom they must not be rotation 
invariant. Thus, they must satisfy the admissibility condition, 

det(Q ai ({R})) + for g a ({R})=0, (3) 

except maybe at singular configurations. We also assume that, except for singular configurations, 
the relation det(7^ b ({J2})) ^ for Q a {{R}) = is fulfilled so that the frame manifold possesses 
a tangent space at {R}. Typically, {R a = 0} is a singular configuration. Furthermore, for the 
two sets of conditions in Q to be compatible the rotational conditions Q a must be translation 
invariant, G a ({ r a + v }) = Ga({ r a}) for any v. This is satisfied by all usual frame conditions (see, 
e.g., In what follows, however, it will be enough to assume only the weaker form, 

N an 

a=l J 

The condition C = 0, on the other hand, is clearly rotation invariant. 

The frame conditions are obviously not unique. Consider the class of reparameterizations, 

3 

^(M) = E^(M)^(M) , (5) 

6=1 

where Vab{{?}) is non-singular on the frame manifold, 

det(V ab ({R})) + for Q C ({R}) = . (6) 
The frame conditions Q' a define the same frame as Q a ctnd we have, 

*3 = E ^rir d ir = E ^v ddl ni idi , Q > ai = . (7) 

iu a n a j ix a j , , 

a=l J J ci,di=0 6=1 

Notice that the reparameterization @ is not necessarily linear in the frame conditions, since the 
coefficients V a b can depend on Q a . For instance V a b = T^-ab satisfies (jSJ by assumption, and leads 
to Q' a orthonormal on the frame manifold, i.e., TZ'^ b = 5 a b- Thus, we could assume without loss 
of generality that frame conditions are orthonormal in this sense. We shall not do so, however, 
because orthonormalizing a set of frame conditions can be inconvenient in practice. 
Another important example is furnished by the reparameterization 

^(W) = E^(W)ft(W). (8) 

6=1 

From © we see that © is satisfied. Clearly, T{ depend non- linearly on the Q a . If Q' a are any 
frame conditions equivalent to the Q a then, in a neighborhood of the frame manifold Q a = 0, they 
can be related by a reparameterization of the form (jSJ). Using ((7J) and (jHJ) we see that the frame 
conditions Ti are invariant under reparameterizations of Q a ^ 

3 3 
mr}) = E ^H{r})G' b ({r}) = E <^c{{r})Gc({r}) = ^({r}) . (9) 

6=1 c=l 
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Other reparameterization-invariant quantities involving the frame conditions Q a can be expressed 
most economically in terms of T{ and their derivatives 



These quantities play an important role in what follows since, as shown below, the frame con- 
ditions Q a enter the Hamiltonian only through T% and their derivatives. This ensures that Ti is 
reparameterization invariant, as all observables should be. 

The previous analysis can be extended to a more general class of reparameterizations of the 
frame conditions Q' a = V a [{Gb}}, involving arbitray functionals V a of Qb which are not singular at 
Qb = 0. Clearly, only the behavior of the new frame conditions Q' a in a small neighborhood of the 
frame manifold Q a = is relevant. In such neighborhood we can expand Va about Q a = 0, thus 
obtaining a reparameterization of the form ©. 



2.2 Body frame transformation and angular momentum 

The transformation relating a configuration {r} of the system in the laboratory frame to the 
corresponding configuration {R} in the body frame defined by the conditions is 

R a = U(r a - r cm ) (11) 

with r cm the center-of-mass position in the lab frame and U = U({8 a }) an orthogonal matrix 
parameterized by three angular variables {#a}a=i- Although our approach and results do not 
depend on any specific parameterization of the rotation group, some parameterization-dependent 
quantities such as the momenta pe a conjugate to 9 a are physically meaningful and play an important 
role in some intermediate calculations. All the information we will need about the parameterization 
of U is encoded in the matrices A and A defined by 

8T7 dU 

C/t = A ai Ti, C/t =\ ai T i , a = 1,2, 3, (12) 

where is the tranpose of U and Tj are the standard generators of the so(3) algebra, (Tj)^ = 
Eijk- The three matrices dU /dO a U\ a = 1,2, 3, must be a basis of so(3) for all values of {0^} if the 
parameterization is to be well defined. Thus, the matrix A a i is invertible and, analogously, so is X a i- 
From the unimodularity of U it follows that U^TiU = UijTj and then, from (|12|) . X a j = A a iUij. 

The frame conditions (^Q) determine the time dependence of {0 a } so that, given a trajectory 
{r a (t)} of the system in the lab frame, we have G a ({Ra(t)}) = G a ({U ({9 a {t)})(r a (t) — r cm )}) = 
for all t. From (|11|) we then have, with M the total mass of the system, 

dRpi /„ m n \ TT dU, 



Substituting (fT3)l into the relation dQ a /dR a j = 0, and using the definition (j2J) for Q and the 
antisymmetry of (dUik/dr a jUik) in i and I, we obtain the relation 

Tj U[k — ^ £jlm Qma q r> ^ n 3 — ^-iim-F manUnj i (1^) 
UT a j Oii an 
J a=l 
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which expresses dU / dr a jU^ in terms of and {Ob}- This expression characterizes the depen- 

dence of U on {r a }, and will be important below, especially in the discussion of angular momentum. 
Some further consequences of (|14j) are discussed in appendix [XJ Together, Q13J) and (|14l) lead to 



dR 



dr aj 



<W JT) Uij + Sikl^lamUmjR/3k , (15) 



which will also be useful below. 

From (|TT|) and dR a i/d9 a = we get 

dr ai dU ki 

-oe; = ^ Uk ^- rc ^ ) - (16) 

We assume that interactions among particles do not depend on their velocities. Classically, the 
momenta pg a conjugate to 6 a is then pg a = dC/dO a = ^2a=i(^ r m/ '90 a )(dC/ 'df a i) where C is the 
classical Lagrangian in the lab frame. Thus, taking 1)12(1 and ()16|) into account, we have 

Pe a = -Kiih ~ Icmi) = ~KjLj with l cra = Mr cm A r cm , L = U(l-l cm ). (17) 

Z cm is the center-of-mass angular momentum in the lab frame, and L the total angular momentum 
about the center of mass in the moving frame. 

In the lab frame the angular momentum operator I satisfies the usual commutator algebra. 
Using dJ, the definition Q of Q, and the unimodularity of U, we obtain, 

N ^ 

^ % dr am 

a=l 

Using (|18|) and the definition (|17|) of L, its commutators can now be computed 

[Li,Ujk] = —i^ijnUnk , [li,Lj] = , [L{,Lj] = —iEijkLk ■ (19) 

The commutators among components of I and L are the same as those for a rigid body, with I the 
angular momentum in the laboratory and L in the body frame. We notice also that [li, R a j] = = 
[Li,R a j] as expected, since in the body frame rotations act only on the angles {9 a }. 

2.3 Internal coordinates 

In order to describe the dynamics we introduce a set of 3iV — 6 internal coordinates {t a }^]~ 6 
defined locally as independent rotation- and translation-invariant functions of configuration space 
t a = ta({ r })- Some consequences of the Euclidean invariance of t a which will be used below are 
t a ({r}) = t a ({R}) and, introducing notations analogous to the first of 0, 

dt dt dt d 2 t d 2 t d 2 t 

^ a ^d^J {R}) = u d^ i{r}h M-ik- 55 d^- {{R}) = ^- (M) (20) 

which can be derived using (jll|) and ()14j) . The body frame, specified by the frame conditions Q), 
fixes relations of the form R a = R a ({t a }) such that the conditions Q are satisfied identically when 
evaluated on R a ({t a }). Thus, the functions R a ({t a }) are a parametric solution to the conditions 
(PQ). Through the relation inverse to (fTT|) . r a = ({Oa^Radtf,}) + r cm , the internal coordinates 
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{tb} together with the orientational and translational ones, {9 a } and r cm resp., give a set of 3N 
local coordinates in configurations space. 

Usually, internal coordinates are given as 3iV— 6 independent functions {t a } of as many indepen- 
dent Euclidean invariants chosen out of the set of all dot and triple products among {r a — fp}ap=i- 
One way of introducing local coordinates in configuration space is to start with a parametric solu- 
tion {R({Q})} to the frame conditions, where {Q} = {Q a }a=i 6 * s a se ^ °f independent parameters. 
These parametric relations can be inverted to give 3iV — 6 local coordinates Q a ({R}) on the body 
frame. Such "inverse" is clearly not unique, however, since at each point {R} in the body frame 
the functions 

3 

Q'a(M) = Qa({r}) + £ Xa b ({r})G b ({r}) + \ a ({r}) • C({r}) 

6=1 

with X a b and A a arbitrary coefficient functions, take on the same values as Q a ({r}). That ambiguity 
can be fixed by imposing additional conditions such as Yla=i 9Qa/dR a = and either 

dQa dQ b 



y- 

/ — ' m a dR a dR a 

a=l 



1 < a < 3 , 1 < 6 < 3iV - 6 (21a) 



or 

J2 RaA 7nr = ' i<fo<3iv-6. (21b) 

a=l 

The first set of conditions fixes the ambiguity because detfK 2 ({#})) / 0, and the second because 
det(Qai({R})) 7^ 0, by assumption. Once the ambiguity has been fixed, to each configuration {R} 
in the body frame there corresponds one and (except at singular points) only one set of parameters 
Q a ({R}). We call those parameters body-frame coordinates. Rotation- and translation-invariant 
coordinates in configuration space can be obtained locally by extending the body-frame coordinates 
by rotation and translation, t a ({r}) = Q a ({R({r})}), with R a ({r}) given by (|TT|) . These t a are 
invariant under Euclidean motions because {R({r})} are. 

The body frame coordinates Q a , considered as local functions of configuration space Q a {{r}) 
are not, in general, rotation or translation invariant. For instance, when the frame conditions are 
linear the Q a ({ r }) can be chosen to be linear functions of {r}, therefore not rotation invariant, 
as in the case of normal coordinates in the Eckart frame. A simple example is given in section 
16.11 The derivation of the Hamiltonian operator and Hilbert-space inner product in terms of linear 
body-frame coordinates and their conjugate momenta has been discussed in detail in p]E|- I n 
this paper we confine ourselves to a dynamical description in terms of Euclidean-invariant internal 
coordinates {t a }. 

2.4 Hamiltonian operator in configuration space 

The kinetic-energy operator in the lab frame is given by the familiar expression 

N 

m a dr a jdr, 



In this equation and in what follows we use units such that h = 1 unless otherwise stated. In terms 
of mass-weighted position vectors r' a = ^/m a r a , fC takes the form of (—1/2 times) a Laplacian 
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operator in Cartesian coordinates in 3iV dimensional space. Given a set {q} = {q\, . . . , <73jv} of 
curvilinear coordinates in configuration space we can write K, as a Laplacian in either one of two 
commonly used forms. Application of the chain rule to (|22|) leads to the standard ordering, with 
all derivative operators to the right of coefficient functions, 



3N „ „ , 3W 



a, 0=1 6=1 

, -l _ 1 dq a dq b - m dr aj dr aj u _ 1 <9 2 <? 



(23) 



— m a dr aj - dr aj - ^ d<? a c% m a dr Qj dr aj 



This form for the kinetic energy has been widely used in molecular physics and leads to expressions 
which are usually simpler than in Weyl ordering. An equivalent expression is 

1 3/V r, r, 

a, 6=1 

Despite their equivalence, (|23|) and 1)24(1 lead to considerably different forms for the kinetic energy 
operator for many-body systems, especially after the momenta conjugate to orientational variables 
are eliminated in favor of the total angular momentum. Those forms, and the relations among 
them, are discussed below in sections |3] and 03 



3 Hamiltonian operator in standard ordering 

The kinetic energy can be expressed in terms of internal coordinates within a given frame convention 
Q by taking {q} in (|23|) to be the union of the set {i a }j^~ 6 with the orientational coordinates 
{@b}t>=i an d the lab frame center-of-mass vector r cm . The kinetic energy term depending on 
r cm trivially decouples from the other degrees of freedom so we simply ignore it in what follows. 
We introduce the notations p a = —id/dt a (a = l,...,3iV — 6), pg b = —id/dOf, (6 = 1,2,3) and 
K, = K,q+K,qi+K,02-, where /Co, the vibrational kinetic energy, does not contain pg a , ICgi contains the 
term quadratic in pg a and those terms linear in pg a whose coefficients involve only first derivatives 
of 9 a with respect to {r}, and KLq2 gathers the remaining terms linear in pg a , with coefficients given 
by second derivatives of 9 a . 

An expression for /Co can be immediately obtained from (|23|) as, 

^ 3JV-6 ^ 3/V-6 

^0=2^ 9abPaPb + 2~ 9bPb ' 

a,b=l 6=1 



dab 



(25) 

E l dt a dt b _ 1 d 2 t b 



ni a 0R a j 0R a j ' m a dR a jdR a j 



where use was made of (|2U|). Notice that g a b and g b depend only on {t a }. 

In order to obtain fCgi we need to specify the dependence of the orientational coordinates {8 a } 
on the lab frame coordinates {r}. Using (J12[) we have 

p. ^ Ulk = ^ —AbmZiml • (26) 

or a j or a j 
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Equating the r.h.s. of iffify to that of ((Til) , we obtain the dependence of 8 a on lab frame coordinates 
in terms of the frame conditions |T|). 

d9 

— — A. ma J~ rnaiUij , 0. — 1,2,3. (^7) 

This expression is independent of the frame-conditions parameterization, as it should be. From 
(|27jl and (|2U|). the corresponding blocks in the matrix fc"^ 1 defined in ()23|) are 

, -i _ 1 dt a d9 b _ ^ 1 gtg ! 

" 1 (28) 



* 1 <96» a <96> b 3 ^ 1 

fc 0ak = E ^~fh^~fh~ = ^rnl^JaKb , ^i/ = E ^cdQ-icQjd = E '^T^ak^ 
a=l m ° otj or aj c d=i a=i m a 

The remaining block fc^ is obtained from k~g by symmetry. The coefficients (j28j) fix the form of 
K,q\. Furthermore, we can eliminate p$ a in favor of the body- frame angular momentum L with the 
aid of (|17|). Taking account of the ordering of operators we get, 

3N-6 3 r, 

o=l - o,6=l " '' (29) 

N 



dt a 

m, „, 



£ — J rn 



m a dRal 



The coefficients of the first two terms in Kq\ depend only on {t a }. There still is dependence on 
{9 b } in the last term, which will cancel against an analogous term in IC$2- We turn to the latter 
next. 

The operator is linear in pg a , with coefficients kg a given by (|2*3*|) with q b substituted by 9 a . 
The second derivative of 9 a is obtained by differentiating both sides of (f2*7|) with respect to r a j. 
Using 

dKrta = dO d dA^ a 
dr ocj £^ dr a j d9 d 

with d9d/dr a j given by (|27|i. and the expression for dUij/dr a j from (|14|>. we get, 



N „ Ar 

« _ \ " 1 OJ ~lai TT _ \ ^ 1 -r ,r- 

^-^ m a or a j m a 



(30) 



Since ^>ai by its definition (|1(J|> depends on {t a } only, using the rotation invariance of t a we can 
write, 

$1 — ^ ^ I n p I" ^inm-F'i man-Flai ] j (31a) 

a=l v ' 
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with 



or- = -^rr(W) = ' (31b) 

" J i-ai ul en ut/ a ulx -ai 

a=l 

(In appendix we rewrite ()31b|) in a completely different way in terms of U.) Therefore, 

a=l a=l 



Thus, /Co, /Cei and /C#2 add up to the total kinetic energy operator, 

3AT-6 

P I. !>.. 1 . 1 4- ^ .. 



' 1 

JC = /C + /Ci + /C 2 , /Ci = y] V ak p a L k + -B k L k , /C 2 = -M^LiLj , (33) 



o=l 



where /Co is given in (|25|). P a fc and are given in ()29|) and l|31aj) resp., and A/^" 1 in $2$. The 
notation used in P3*|) is such that /C n , n = 0, 1,2, depends on the n th power of the body-frame 
angular momentum. The coefficient functions in /C depend on {t a } only, as expected by rotation 
invariance. Clearly, /C is invariant under frame-conditions reparameterizations, and independent 
of the parameterization of the rotation group used to define the orientational variables {9 a }- The 
second term in IC± is purely quantum mechanical, since B k oc H as can easily be checked by 
dimensional analysis. The origin of this term lies in the operator ordering, in the same way as 
other orderings (such as e.g., Weyl ordering) give rise to quantum potentials pQ. The form of B k is 
further simplified by the fact that the first term in the parentheses in (j31a|) vanishes for the most 
usual choices of frame such as the iV-body Eckart frame, and also in the examples in sections El 
and below. 



4 Inner product. Matrix and rigid-rotator Hamiltonians 



Due to translation invariance the lab-frame wave function can be factored as ^(o)({ r }) = V^o) ({ r » — 
r cm }) exp(ifc cm -r cm ), with the subindex (0) indicating lab frame. Starting with the canonical inner 
product in the lab frame and changing variables to {t a }, {9b}, r cm we get, 



3N-6 



(34) 



0=1 



6=1 



Here we already integrated over Ton > obtaining a momentum-conservation 5 function which, we 
omit. The Jacobian J can be expressed in terms of internal coordinates by means of the relation 



j 



det 



dq a 
dr, 



-1 



det 



dq a 
dR 



-i 



(35) 
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where |A| = det(A a j), with A a j defined in (jTl)) . and dq a /dR a j = Ukjdq a /dr a k and we used (|27|l . 
1/3 in ()35|) is the absolute value of the determinant of the 3 N x 3N matrix 



/ dti/dR lx dh/dR 1Y 

dt 3N - 6 /dR lx dtsNs/dRiY 
—Fwx —F\\y 

—F31X —F31Y 
dr cm x/dr lx dr cmX /dri Y 



• • dh/dR NZ \ 

.. dt 3N _ 6 /dR NZ 
—Finz 

—F3NZ 
■ ■ dr cmX /dr NZ 



(36) 



\ dr cmZ /dr lx dr cmZ /dr 1Y ... dr cmZ /dr NZ J 

The minus signs in the three middle rows are of course unimportant in ()35|) . Notice that the 
quantities Ti a j are completely determined by the frame conditions, and are usually much simpler 
in form than gradients of Euler angles. The Jacobian J as given in (|24jl . which is proportional to 
3 , can be computed with the derivatives (|27|) in terms of the matrix g a i from (|25|) . The procedure 
in this case is closely analogous to the case of linear frame conditions discussed in jTj . We will not 
dwell on that calculation, whose result and its derivation have been considered in [3]. With our 
notation we have, 

J=Urnl /2 J = M^\A\ 1 -^, (37) 

a=l 1^1 

where M = Y2 a m °' M is the body-frame inertia tensor and | M\ its determinant and \g\ = det g a b 
with g a b from Below we denote dVg = n&=i d9 u \A\ the invariant measure on SO(3), with total 
volume Vq = 8n 2 . 

Since K, commutes with (Z — Z C m) 2 = L 2 and (l z — l C mz) we can choose ^(0) = V'(o)^m({ r a ~~ r cm}) 
to be an eigenfunction of those operators. The body-frame wave functions are then 



iP im ({R}, {0 a }) = i> mm ({U\{e a })R}) = V2£ + l £ iP ms ({R})D e * s ({O a }) (38a) 



s=-e 



with 



Dt m {{O a }) = j d 2 ey4,(e)y im ([/({^})e) 



(38b) 



where e is a unit vector varying over the unit sphere and D^ n , m ({9 a }) the irreducible matrix 
representing the rotation U({9 a }). In terms of the wave functions we have, 



i ,. 37V-6 



<*/J*W> = Smm'Ve J II dt aJr {0 )J{R})^0), n ({R}) ■ 



(39) 



n=-i a=l 



The action of the lab frame angular momentum on the body-frame wave functions (|38j) is given 

by, 



s,q=-l 



qui 
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where C\ is the standard angular momentum Hermitian matrix in the representation of irreducible 
tensors of order £, [S] 



C ^\m = l d ^ Y ll^)^p\^rn{e) = £ [KaiM^^j ^ £>L(K» • (41) 



The matrix A^daa}) in ()41l) is as defined in (|12|) . Analogously, the body-frame angular momen- 
tum operator acts as, 



L^uiRhiOa}) = -VWTT Yl AoMi R }) (cf ] ) ar Di: q {{0a}) ■ (42) 

s,q=-£ 



sq 



Thus, we can represent K, by means of its matrix elements between angular-momentum eigenfunc- 

((■) 

tions in terms of the matrices C\ . We need consider only matrix elements between wave functions 
with different "radial" quantum numbers, but the same angular dependence, so that, 

i 

— [dV e $ m ({R},{6 a })}Ci Hrn ({R},{9 b })= <T r mp ({R})£ pq i> m ({R}) 



Vg 

£ vq - JC 8 m - ^ V akVa (Cf) - l -B k (Cf) + 1 -NJ {C?) (C 



P ,i=-i (43) 

3N- a 



a=l 



In this equation the operator /Co and the coefficient functions V> ak , B k and M i • 1 are as in (|33j) . 
Instead of the body frame wave functions (|38|) we can introduce the alternative basis, 

M{R},$) = E ^(o)im({Rm* m (e) ■ (44) 
m=-l 

These wave functions, not eigenfunctions of l z , depend on a unit vector e representing a fictitious 
rigid rotator of total angular momentum I. In terms of <j>£ we have, 

y J dV e fciW, {da})L p ^ ls ({R}, {9 b }) = 5 rs J d 2 e^({R}, e)S p &{{R\, e) , (45a) 

1 d 

Sp = -ZpqreqTT- ■ (45b) 



i w q de 



r 



The factor 5 rs on the r.h.s. of <j45a|) is due to the fact that L p commutes with l z (see (|19|l). 
Expression ()45b|) for S p is not affected by the constraint e • e = 1, that is, the derivatives can be 
computed without taking that constraint into account, as can be easily checked. Alternatively, the 
operator S p can be expressed in terms of the spherical angles 9, <f> of e and derivatives with respect 
to them. Therefore, from (|33j) and Q45aJI we have 

y f dV e $ m ({R},{Oa})Wlm({R},{0b}) = I d 2 e^{{R},e)KU{R}^) (46) 

with K, having the same form as K, in (jSSJ), but with the angular momentum L replaced by S 
as given by (|45b|) . Actually, the product LjLj in (|3*3*|) is mapped into SjSi, but that product is 
contracted with A/^" 1 which is symmetric. This rigid-rotator formalism, based on the wave functions 
Q44|) and the Hamiltonian (|46[). is a useful alternative to the matrix formalism based on (|38|1. (|41jl 
and (|43|) . It appears naturally in the gauge- invariant approach of pQ. 
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5 Other forms for the Hamiltonian 



Other forms for the many-body Hamiltonian in a body-fixed frame, based on expression (|24[) for 
the Laplacian, have been given in the literature (see, e.g., [3J|!J] and references therein). In this 
section we discuss the derivation of the kinetic energy operator in the form Q24JI from the point of 
view of frame conditions and establish relations among these results and those of section Our 
notation follows that of P|. 

In order to express the kinetic energy in the form (|24[) it is convenient to write the matrix k~ h l 
in a form different from that used in sec. El Defining, 



8Rni dR n 



lab — 

a=l 



and using the chain rule we get, 



FH- 3N ~ 6 N f)F! ap 

6=1 p=l 

Substituting (|4*8|l in the definition (|25|) of g a b, and using the derivatives (|15jl. we get, 

fafe = ^ + X] hacKd^djK^ 1 ( 49 ) 
c,d=l 



with A/"j • 1 defined in (j2BJ) and 

N dR 
a d = Y^ mpRp A — £ . (50) 

/3=1 a 

Several relations between g a0 and h a b, and between M" and A/", analogous to (jlU)) are summarized 
in appendix U3 Eq. (J3HJ) fixes the form of k~ b x in (J2U for 1 < a, 6 < 3iV - 6. 

Similarly, we can obtain a compact expression for the off-diagonal blocks of k~. . Using (|48jl 
and (|15|) together with the frame conditions Q a and their translation invariance, we obtain, 

3JV-6 

Ki = £ Kt**Nj*x& > (51) 

c=l 



to be compared with the corresponding expression in (|28[). The block fc^g is as given in (|28[1. As 
in sect.|3]we omit here for brevity the terms involving the center-of-mass degrees of freedom, which 
are dynamically trivial. With these expressions for k~^ we can compute its determinant 1/J 2 , by 
factoring the matrix appropriately. We again omit the details and state the result 

J = M 3 / 2 |A||AA| 1/2 |/i| 1/2 , (52) 

with \h\ = det(h ab ) and \Af\ = det(AA). Comparing Q% with (J27J) yields \M\/\g\ = \h\\Af\ 0. 

The matrix k~^ and J are all we need in order to obtain the kinetic-energy operator K, from 
(|24l) . We can, however, eliminate all dependence on orientational degrees of freedom by means of 
the well-known relations (see and refs. therein) 

^XX^mE^A-IAI. (53) 

a=l a 1 1 6=1 
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With this, we finally get, 



1 37V-6 

(54) 

, 3JV-6 \ / 37V-6 x 

+ WWW 2 [ Li ~ ^ Phh ™ avi ) \ Af \ 1/2 \ h \ 1/2 ^ 1 l L i- E h ~>w 

Notice the ordering of operators in ()54|) . K, can be Weyl ordered most easily after performing 
a transformation of the form J7C 1/J, leading to a quantum potential term. Weyl ordering is 
considered in detail in the case of linear frame conditions in pQ, and the same procedure can be 
applied to the case of general frame conditions discussed in this paper. The expression for the 
quantum potential, however, seems to us to be too complicated to be useful in practice so we 
omit the results. Other orderings are of course possible, although with similar caveats about the 
associated quantum potentials. That is an advantage, from our point of view, of the standard 
ordering given in section 12 

We can rewrite /C as given by (f51|) in terms of g a \, and Atf" -1 , instead of and A/" -1 , by using 
relations (JUJ), (JHZJ) and dS2J, and (|B~6|) . to find 

1 1 lol 1 / 2 3N ~ 6 \M\ 1 / 2 

K = o M ^ L i + 9^1/2 E (P° " A ^ TTm9ab (Pb ~ A bj L 3 ) . (55) 

Here we defined [3] 

A a = M~ l a a , 1 < a < 3N - 6 . (56) 

In the form ((53|) all dependence of /C on frame conditions is implicit in the relations R a = R a ({t a }), 
which enters K, through g ab and A a and also through M^ 1 when expressed in terms of internal 
coordinates. It is interesting to point out that for the most commonly used frames the 3x3 matrix 
Q a i defined in © is much simpler to invert than M, and that is the only matrix inversion needed 
to obtain K, as given in (|3*3*)) . 

Comparing the expressions (|54j) and (|55|) with ()33|) we can obtain relations among their coeffi- 
cients. The equivalence of the purely vibrational terms in (|33|) and (|54j) is immediate once we take 
into account ()B.8|) and the equivalence between the two standard forms for the Laplacian ()23|) and 
(|24|) . Similarly, the terms quadratic in L in (|33|) and (|54|) are obviously equal and equivalent to 
that in (|5o]) by (|B.6|) . Notice that ftf^j , which is usually defined 3 as in ()B.4|) or ()B.6|) . can be 
compactly expressed in terms of frame conditions by our definition (|28|) . 

Equating the terms linear in L in (|3*3~|) . and (|5oT) we get the relations, 

3iV-6 3JV-6 
V dq = -M' 1 ^ gdd'dd'r = -Nqr 1 ^ h dd' a d'r (57a) 



d'=l d'=l 



I 



1 | 5 |l/2 / |M|1/2 



3N- 



with and ,6 g defined in (|2"§)) and (|30f> . respectively. These relations can be proved directly, 
providing a consistency check on our results. An important consequence of ()57a|) is that it allows 
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us to write A a , at least locally, in terms of frame conditions, 

3JV-6 3^-6 

A bi = - ^2 dbc^ci , or a bi = -Nik ^2 h ^ v a , a = 1, . . . , 3N - 6 , (58) 

c=l c=l 

with T> c i given by ()29j) . Notice that these relations cannot be obtained from (|49|) or the equalities 
in appendix [BJ which always involve a a or A a quadratically. Through (|58jl we can write any 
expression involving the gauge fields A a (3] in terms of internal coordinates and frame conditions. 

6 The case N=3 

We consider here the case N = 3 both as an example and a verification of the foregoing, obtaining 
the Hamiltonian in two different body frames, one defined by linear conditions and the other by 
quadratic ones. We choose internal coordinates t\ = pi, t 2 = P2 and t 3 = 9 which are standard in 
molecular physics, 

Pi = \ri-r 3 \, P2 = \r2~r 3 \, cos 6 = — *— (n - r 3 ) • (r 2 - r 3 ) . (59) 

P1P2 

with conjugate momenta denoted by p a , a = 1,2,3. We define also the reduced masses l/mi3 = 
1/mi + 1/7713 and analogously 711,23. 

6.1 Linear frame conditions 

A linear body frame with origin at the center of mass can be defined by choosing the Y axis 
orthogonal to the plane of the system, Y oc (r 2 — r%) A (7*1 — r 3 ), the Z axis along r\ — r 3 , and 
X = Y A Z. The frame conditions are then, 

N 

C = ~^ ra = ' £1 = r iz - r 3:r = , £2 = rij, - C y = , £3 = r 2y - C y = . (60) 



£*=! 



Notice that £/ a are written so they are explicitly translation invariant. From (|25| ) and (|59|) we get, 

1 cos sin 6 

9n = j 512 = , 913 = 

mi3 m 3 m 3 p2 

1 sine 1 1 2 cos 6 

922 = , 523 = , 533 = 2 H 2 (61) 

71123 m-iPi rn 13 pf m 2 3P 2 ™3PiP2 

2 2 . / 1 1 

5i = , 52 = , 53 = cot 6M j H ^ , 

mi3/9i "123/92 V^isPi m23PiJ m 3 pip 2 sm 

and the Jacobian entering the inner product l)39|) is found to be J = p^p^smO. The coefficients 
(j61|) fix the form of /Co as given in l)23|) . With the frame conditions (|60j) and the internal coordinates 
from (jTUJ) we obtain, 

= — y , = o , :Fi3 = — y , 

pi pi 

T 2X = —X , ^22 = , JF 23 = —X , (62) 

pi Pi 

F 31 = _^i> jr__ = 1 Y ^33 = _P1-^CQS^ 

pi p 2 smt) pip 2 smt) 
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with the notation tF aa = (J~'aai,3 r aa2,3 r aa3)- We omit the details of the calculation ofjsf^j 1 , V a k and 
Bi (see (|28|l. (|2*^|) and (|30j). resp.). Rather, we give the result for /C, from which those coefficients 
can be read off. The kinetic-energy operator for this system is given by (|33|) as K, = ICq + /Ci + K2, 
with /Co resulting from (|61|) and with 

_ sin# / 1 cos# \ i / 1 cos# \ 

M = p 2 ^y + 2 + . fl 2 L ^ 

"i3Pi V m i3/0i rn 3 p 1 p 2 J 2 sin \m 3 p 1 p 2 m V3 p{) 

^2 = ^-^L 2 x + U J——) {L X ,L Z } + _L^L!- (63) 

2mi 3 pf 2 \m3p1p2smfc' mispf tan 9 J 2mi 3 p\ 

+ 1 (( 1 + _i L_ _ 2cosg ^ L 2 

2 V\ m i3Pi m23/)|ysin 2 6l mi 3 /3f m 3 pip 2 sm 2 9j Z ' 

These results agree exactly with those of [2] once we take into account that the kinetic operator 
defined there is p\p 2 K. l/{p\P2) in our notation. 

In this example, since the frame conditions are linear, we can choose a set of linear body-frame 
coordinates satisfying (|21aj) . We set Qi = R\z — R3Z, Q2 = R2Z — R3Z, Q3 = R2X — R3X- 
These coordinates Q a can be extended to all of configuration space by linearity, yielding a set of 
non-rotation-invariant coordinates Qi({i*}) = T\ z — r 3z , etc. In order to extend them to rotation- 
invariant internal coordinates we express them in terms of scalar products of body-frame position- 
vectors. Such procedure leads to a set of coordinates equivalent to (|5l3|) . which with the same 
notation are written as pi, p 2 cos9, p 2 sm9. 

6.2 Quadratic frame conditions 

Another frame for the three-body system used in the molecular-physics literature is defined as a 
modification of the previous one, choosing the Z axis to bisect the angle 9 between n — r 3 and 
r 2 — r 3- The frame conditions are as in 1)60(1 . except that Q\ now takes the form 

Si({r}) = (ri x - r 3x )(r 2z - r 3z ) + (r lz - r 3z ){r 2x - r 3x ) = . (64) 

This frame differs from the body- frame of section UTTl bv a time-dependent rotation in an angle 9/2 
around the Y axis. Since internal coordinates are rotation invariant, g a b and g^, and therefore also 
/Co, are as in (jfilj) . The Jacobian J also remains the same as above. 
The modified frame conditions lead to, 

F\\ = -- — " — Y > ^21 = ^- cos °-X-^- sin , T 31 = 1 Y 
2picos| 2pi 2 2pi 2 2pism| 

Fi2 = - ~ 1 e Y , ^22 = ^cos^X + -^sin^Z, ^32 = - 1 e Y (65) 
2p 2 cos| 2p 2 2 2p 2 2 2/9 2 sm I 

^13 = -- * . -^23 = -^ COS-X-- Sm-Z , ^33 = - — g * , 

2P1/J2COSI 2p + 2 2p_ 2 2pip 2 sm| 

where the notation is as in Q62(l and l/p± = il/pi + l/p2- The kinetic energy is then fC = 
/Co + /Ci + /C2, with /Co given by (|HT|) and with 

sin0 / 1 1 \ 1/1 1 \ /cot0 \ _ 

£1 = q —Pi P2 + - 2 2 "^F - + P 3 L ^ 

2m 3 VP2 Pi J 2 \m 13 pf m 23 p z 2 ) \ 2% ) 
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8 cos 2 1 \mi 3 p( ' m 2?> p 2 ' rn 3 p 1 p 2 J X 8 \m 13 pl m 23 p\ 



^ = 7—nf —2 + z—2 + z— ) l2 x + o\ —2 + z—s + z— L l (66) 



I f I 1 2 \ r2 1/1 1 \ r , 

+ ~ . o 6> 2 2 L 2 + A ■ a 2 2 X^Xi^Zi ■ 

8sm 2 | \ m i3/0i rn 23 p^ ™>?,P\P2 ) 4sm0 V"W2 m nPiJ 

This expression for /C agrees with the result given in as corrected in j2J 3 taking into account 
that their operator corresponds to pip 2 IC \/{pip 2 ) in our notation. 

7 The case N=4 

As a further example we consider in this section a four-particle system. Our choices of frame 
and internal coordinates below are appropriate for a system with the topology of the formalde- 
hyde molecule, though the results are also applicable to other systems for which those choices are 
not singular at the equilibrium configuration. The vibrational Hamiltonian for the formaldehyde 
molecule has been given, in the Born-Oppenheimer approximation, e.g., in [2] (see section 4.3 and 
appendix A.) Other explicit results for four-body systems are given in 0] and references therein. 
Here we use a set of internal coordinates which, combined with the general results given above, 
greatly simplify calculations and lead to moderately simple results for the total Hamiltonian, in- 
cluding rotation and vibration-rotation terms. One drawback of our coordinate choice, however, is 
that it also results in a complicated expression for the inner product. This section is not meant as 
an exhaustive kinematic analysis of the four-body problem, but rather as an example of the results 
given above. 

In this section we label the particles with capital letters, a = A, . . . , D. In the case of the 
formaldehyde molecule D would refer to the carbon atom, A to the oxygen, and B and C to the 
hydrogen atoms. We choose a frame with origin at the center of mass whose Z axis lies along 
Rad = Ra — Rd, and the Y axis is defined by the condition that Rcd lies on the Y Z coordinate 
plane. This choice of frame is singular when Rcd is parallel to Rad- The rotational frame 
conditions are, 

Qi = tadx = , Q 2 = r A Dy = , Q 3 = r C Dx = . (67) 

The frame is completely determined by (|67j) together with the auxiliary conditions Radz > and 
Rcdy > defining the direction of the axes. From (jBTjl and © we get, 

/ p-^ D ii- D x \ I Radz \ 

Kb=\ Pad > Qai = -Radz (68) 
V p D l Pcd) \ Rcdz -Rcdy J 

with 1/ Pad = 1/ PA + 1/ Pd and similarly for the other reduced masses. The frame conditions are 
therefore singular when det(Q) = —R 2 adz Rcdy = 0, i.e., when Rad and Rcd are parallel, or 
either one vanishes. With the matrix Q in (|68|) from (|10|) we obtain 

^1A2 = — t; i FlD2 = T, ) F2A1 = T, ' F2D\ 



Radz Radz Radz Radz /gg\ 

Rcdz ^ 1 ^ Radz - Rcdz 

•^3A1 — -5 5 , -T3C1 — -"5 , -^3£>l — —5 5 , 

tiADZtiCDY tiCDY tiADZtiCDY 
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all other Ti a j vanishing. In turn this leads to 

Ki= — A/23 1 = — ^ =M£, 

HadRadz hadRadz^cdy HdRadzRcdy 

1 KT-l R CDZ 1 n -RCDZ 



(70) 



AT- 1 - - AT 1 

■/von „c, , JVqi 



HadR\dz ' "" ^adR\ dz R 2 cdy VcdRcdy HdRadzRcdy 



+ 



and the remaining components vanishing. With AA _1 from 1)701) . the rotational kinetic energy /C2 
(f33|) is completely determined. 

Our choice of internal coordinates is motivated by calculational simplicity. We introduce 
translation- and rotation-invariant internal coordinates depending polynomially on the position 
vectors, 

h = r AD , to = r 2 CD , t 3 = r CD ■ r A D , ^ 
U = r B D ■ rAD , h = r BD ■ r C D , h = r BD • r C D A r A D ■ 

1/2 1/2 

The ranges of variation for the first three coordinates are t\ y2 > 0, \t 3 \ < t 1 t 2 , whereas the 
last three can take any real value. In the frame defined by conditions 1)67)1 and the associated 
suplementary conditions, relative particle positions are given by 

R AD = VhZ , R CD = ^|EI y + J* z , 

V*i V'l ( 7 2) 

Particle position vectors Ra,b,c,d can of course be found from (|72j) together with the center-of-mass 
condition. 

The vibrational kinetic energy JCq in standard order is determined by the coefficients g a b and 
gb in 1)25)1 . Due to the polynomial nature of i a , the results for 55 are remarkably simple 

6 6 6 

5i = , 92 = , 93 = 94 = 95 = , 96 = . (73) 

HAD HCD HD 

The expressions for g a b are unavoidably more complicated, even though their dependence on R a 
is polynomial. Expressing g a b in terms of internal coordinates t a we get, taking into account its 
symmetry, 

4 4 2 2 2 2 2 
5n = h , 912 = — h , gi3 = 1 3 H h , gi4 = U H *i , 915 = — (*3 + ^) , 

MAD M-D MAD HD HAD HD HD 

2 4 2 2 2 2 2 

516 = h , 922 = 1 2 , 923 = h H *2 , 924 = (*3 + h) , 925 = *5 H *2 , 

MA MCD MCD MD MD MCD MD 

2 112 11, 

926 = — *6 , 933 = h H ti H 1 3 , 934 = 1 5 H (ti + 1 3 + U) , 

Mc had hcd hd had hd 

11 112 

935 = U H (t 2 + *3 + h) , 936 = , 944 = R B d H *i H *4 , (74) 

Hcd Hd Had Hbd Hd 

1 1 2 1 1 2 2 

945 — *3 H \R-BD + *4 + *5j > 946 — , 955 — 12 H "BD "I *5 j 

MBD MD Hbd Hcd Hd 

956 = , 966 = — (R BD t 2 - tf) H (*2*1 - *i) "I (R-BDh - *1) • 

MA Ms Mc 
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Here we have left R\d indicated for convenience, its expression in terms of internal coordinates is 
given by ((721) . 

The vibrational-rotational coupling term K\ is given in ()33|) in terms of the coefficients P a fc 
and Bk (see (J2§)) and ()31a(l ). The expression for V a k can be written most compactly in terms of 
position vectors. Its non-vanishing components are, 

2 RcDY 1 d -r. 1 R BDY ^ 1 R BDX 

L>2\ = 75 , ^31 = RCDY , LJA1 = 75 , LJA2 = 75 , 

HD tiADZ HAD HAD i^ADZ HAD &ADZ 

P 1 RbdxRcdz 1 Rbdx p 1 Rcdy + Rbdy p _ 1 Rbdx 

Had RadzRcdy hd Rcdy ' hd Radz hd Radz ' 

P _ J_ RbdxRcdz _ 1 Rbdx ^ __ J_ r bdx r cdz 
Hd RadzRcdy hcd Rcdy ' /M -Radz 

r> - 1 /re a R \ Ti - r cdz ^ , 1 Radz Rbdy 

L>62 — 75 (-KbD A Hcd)X , -^63 — 75 ^62 H 75 • 

HaRadz Rcdy hc Rcdy 

The coefficients Bk, on the other hand, acquire a very simple form because the first term in (|31a|) 
vanishes, leaving only the contribution from the second term, 

Bl = 152—75 ( — r cdz + —Radz] = . 1 ., ( + —) , (76) 

R ADZ R CDY V HAD HD J a/ tlt 2 ~ t| V HAD h HD J 

and B 2 = = B 3 . 

Finally, the Jacobian J in the inner product (|39|1 can be computed to give, 
L = A(R AD A Rod) 2 + ^{Rad A Rod) ■ (Had a + Rbd a .Rcd + «cd a .Rad) • (77) 

J is singular at .Rad A .Rcd = 0, as expected from our choice of frame and internal coordinates. 
This singularity, together with the somewhat involved integration limits resulting from (|71[) . make 
the expression for the inner product computationally cumbersome. For systems whose equilibrium 
configuration is far from the singularity, however, the contribution from that region should be 
strongly suppressed by the wave functions in (|39j) . 



8 Final remarks 

In this paper we derived the body-frame Hamiltonian for a system of N particles in terms of 
frame conditions and internal coordinates. Obtaining the Hamiltonian in terms of frame conditions 
instead of Euler angles and the inertia tensor and their derivatives leads arguably to computational 
simplifications. All frames used in applications are defined by polynomial conditions, usually of 
first or second degree. The coefficients M, T> and B in the kinetic energy operator (|33|) . given 
by algebraic expressions in terms of first derivatives of those frame conditions, can be efficiently 
evaluated with symbolic computer algorithms or, depending on N and the internal coordinates 
t a being used, even by hand. In particular, there is no need to invert the inertia tensor, or to 
compute its determinant or that of the vibrational kinetic tensor g a b (|25|) ■ Similarly, neither those 
determinants nor derivatives of Euler angles are required for the computation of the volume element 
in the quantum inner product as given by 1)35)1 and 1)36)1 , and the derivatives of internal coordinates 
involved in 1)36)1 are evaluated only at the frame manifold. Those simplifications should be more 
apparent the larger the value of N. Furthermore, given a set of internal coordinates {t a }, it is 
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straightforward to compute the Hamiltonian in different frames by changing the conditions Q a , as 
illustrated in the examples of section H3 

The Hamiltonian is given in standard order in (|33|) and in the alternate forms (|54[) and (|55|) , 
Comparing those three forms leads to some useful relations among their coefficients, in particular 
the expression (|58jl for the gauge field A a in terms of frame conditions. One advantage of the 
standard-order form (|33f) is that it is known to be equivalent to a path-integral formulation in 
phase space with post-point discretization, whereas for the undefined orderings of (|54|) and (|55|) 
the path-integral equivalents are in principle not known and have to be constructed. In section 
HI in connection with the quantum inner product in the body-frame, we discuss two equivalent 
representations for the angular momentum operators and the kinetic energy. Namely, as irre- 
ducible matrices acting on ^-component wave functions, (|43[) . and as differential operators acting 
on rigid-rotator wave functions, ()46|) . The rigid-rotator representation, which can be a convenient 
alternative to the matricial one for some computations, appears naturally in the gauge- invariant 
approach of p. 

In section |2] we discuss frame conditions from the point of view of their admissibility and 
reparametrizations. Not discussed in this paper is the problem of frame singularities. From Q we 
see that the singular points on the frame manifold are determined by the equations det(Q a i({R})) = 
= G a ({R}), which are polynomial in R for polynomial Q a . Such algebraic formulation of the 
problem might be useful in the study of frame singularities for larger values of N. 

Applications of the approach presented here to the analyisis of systems with N > 3 are currently 
in progress, and will be discussed elsewhere. 



References 



[1] 


A 


0. Bouzas and J. Mendez Gamboa, J. Phys. A 37 (2004), 6773. 


[2] 


N 


C. Handy, Mol. Phys. 61 (1987), 207. 


[3] 


R. G. Littlejohn and M. Reinsch, Rev. Mod. Phys. 69 (1997), 213. 


[4] 


A 


V. Meremianin and J. S. Briggs, Phys. Rep. 384 (2003), 121. 


[5] 


C. 


Eckart, Phys. Rev. 47 (1935), 552. 


[6] 


.1. 


D. Louck and H. W. Galbraith, Rev. Mod. Phys. 48 (1976), 69. 


[7] 


.1. 


Mendez Gamboa and A. O. Bouzas, J. Phys. A 36 (2003), 7061. 


[8] 


D 


M. Brink and G. R. Satchler, Angular momentum, Oxford U. Press, New York, 1993 


[9] 


.1. 


K. G. Watson, Mol. Phys. 15 (1968), 479. 


[10] 


S. 


Carter, N. C. Handy, and B. T. Sutcliffe, Mol. Phys. 49 (1983), 745. 



A Remarks on the body-frame transformation 

Equation ((HJ) for dU/dr Q gives a relation between the dependence of £/({#{,}) on r a and the frame 
conditions Q a . Notice that on the r.h.s. of (fTl]) the dependences on {t a } and {6^} are completely 
f actor ized, with Tmcm depending only on {t a }. We can rewrite (fT4"|) as, 

•Fiotj = "^^ikl~^ Ui m Uj n . (A.l) 

2 dr an 

Differentiating both sides of (|A.1|) and using ()2U|) and (|31b|) we obtain, 

dFjgj _ 1 d 2 U kn 

-gjf- ~ 2 e iki Q g Ui n . (A.2) 
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Notice that the l.h.s. of ()A.1|) and (|A.2|) depend only on {t a }. The l.h.s. of (|A.2|) appears in the 
coefficient £>/ defined by (|31aj) . 



B Some useful identities 

In this appendix we gather some useful relations analogous to (|i9*)l 3j. The last of these, (|B.9|) . is 
a closure relation which must hold on the frame manifold except at singular points. 



fJa 



-1 



h a b ~ aaiabjM^ 1 (B.3) 



3^-6 



Mij = Mij - h J a a a dj (B.4) 

c,d=l 

JV-6 3iV-6 

^2 9abau = Mij-^jk 1 h ab a bk (B.5) 



6=1 6=1 
3iV-6 



AAt 1 =M: 1 +^ Ari5«6^- (B.6) 



a,6=l 
3iV-6 



hj = 9ab~ Y 9acA ci NijA dj g db (B.7) 

c,d=l 

p^l = det(h ab )det(M) (B.8) 

"hMi* =^K ab —— + ^ rn,—h cd m p — + * iJfc _ (B.9) 
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